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An effective-mass theory with a deformation-induced (an axial) gauge field is proposed as a 
theoretical framework to study graphene edge. Though the gauge field is singular at edge, it can 
represent the boundary condition and this framework is adopted to solve the scattering problems 
for the zigzag and armchair edges. Furthermore, we solve the scattering problem in the presence of 
a mass term and an electromagnetic field. It is shown that the mass term makes the standing wave 
at the Dirac point avoid the zigzag edge, by which the local density of states disappears, and the 
lowest and first Landau states are special near the zigzag edge. The (chiral) gauge theory framework 
provides a useful description of graphene edge. 



I. INTRODUCTION 

The graphene edge has attracted much attentionjiri^ 
because it is the source of a wide variety of notable phe- 
nomena. For example, the zigzag edge possesses localized 
edge statesiSi"— The edge states enhance the local density 
of states near the Fermi energy As a result, the spins 
of the edge states may be polarized by coulombic inter- 
actioufii Another type of edge, the armchair edge, does 
not support edge states. The zigzag edge is the source of 
intravalley scattering, while the armchair edge gives rise 
to intervalley scattering. The transport properties near 
the armchair edge may differ significantly from that near 
the zigzag edge^i^^ however, the reason for this variety 
is unclear. 

The Schrodinger equation is a differential equation; 
therefore, an appropriate boundary condition should be 
imposed on the equation. The boundary condition is 
sensitive to the situation of the edge, while the local dy- 
namics, as described by the Schrodinger equation, are the 
same everywhere in a graphene sample. The wave func- 
tion and energy spectrum are dependent on the boundary 
condition. In this sense, the boundary condition is the 
origin of the varietyi^^"— In this paper, we attempt to 
construct a theoretical framework in which the edge is 
taken into account as a gauge field, and not as a bound- 
ary condition for the wave function. We show that the 
framework is useful to obtain and understand the stand- 
ing wave and edge states. 

This paper is organized as follows. In SecUll the quali- 
tative features of the reflections from the zigzag and arm- 
chair edges are shown using the kinematics for elastic 
scattering. In Sec. IIIII a general form of the electronic 
Hamiltonian is given for a graphene sheet with edges. In 
Sees. llVl andlVl the scattering problem is solved for both 
the zigzag and armchair edges, and the standing wave so- 
lution is obtained. A discussion and summary are given 
in Sec. IVll 



II. REFLECTION OF PSEUDOSPIN 

In the inset of Fig. [U we consider the zigzag edge 
parallel to the x-axis, by which translational symmetry 
along the y-axis is broken. Thus, the incident state with 
wave vector {kx,ky) is elastically scattered by the zigzag 
edge, and the wave vector of the reflected state becomes 
{kx, —ky). In contrast, the armchair edge parallel to the 
y-axis breaks the translational symmetry along the x- 
axis, so that the wave vector of the reflected state is 
{—kx,ky). The Brillouin zone (BZ) is given by 90° ro- 
tation of the hexagonal lattice, so that for the incident 
state near the K point in Fig.[Tl the zigzag edge reflected 
state is also near the K point, while the armchair edge 
reflected state is near the K' point. Therefore, scattering 
by the zigzag edge is intravalley scattering, while that by 
the armchair edge is intervalley scattering. 

Pseudospin is defined as the expected value of the Pauli 
matrices (Jx,y,z with respect to the two component Bloch 
function. The pseudospin provides information concern- 
ing the relative phase and the relative amplitude between 
the two components of the Bloch function, and it can 
be used to characterize scattering at the edges The 
Bloch function of the conduction state with wave vector 
k = {kx,ky) is given by 
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where /k = J2a e , denotes the complex conjugate 
of /k, and Ra (a = 1, 2, 3) are the vectors pointing to the 
nearest-neighbor B atoms from an A atom [see the inset 
of Fig. [1] . The pseudospin is then given by 
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(a,)k = 0, (2) 
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FIG. 1: The pseudospin vector field in graphene BZ. Note 
that this field is for the conduction band and the pseudospin 
field for the valence band is given by reversing the direction 
of each arrow. The singularities in this pseudospin field cor- 
respond to the K or K' points. Appendix |X] outlines why 
the pseudospins at the (three) equivalent K (K') points are 
not identical, [inset-top] The hexagonal unit cell of graphene 
consists of A (solid circle) and B (open circle) atoms. The 
xy coordinate system is fixed as shown. The vectors ai and 
a2 are primitive translations. The length of each of these is 
a {a = VSucc, where Ucc is the C-C bond length), [inset- 
middle] The vectors bi and b2 are reciprocal lattice vectors 
defined by ■ hj = 2nSij . [inset-bottom] The vectors are 
expressed as Ri = accBy, R2 = — (\/3/2)acce2, — (l/2)accey, 
and R3 = (\/3/2)acce2, — (l/2)accej;, where e^; (e^) is the 
dimensionless unit vector for the a;-axis (y-axis). 



where 



(3) 



The pseudospin, {{ax)-k, {<7y)k, (o'z)k), may be regarded 
as a two-dimensional vector field, because {az)k = 0. The 
arrows in Fig.[T]show the pseudospin field, {{(Jx)k, {o'y)k)- 
{cry)k is proportional to Ini[/k]; therefore, the angle of 
each arrow with respect to the fc^-axis represents the rel- 
ative phase of the Bloch function between A and B atoms. 
For example, at the F point k = in Fig. [TJ the arrow is 
pointing toward the negative fc^^-axis. This implies that 
the wave function forms an antisymmetric combination 
with respect to the A and B atoms, which can be checked 
by setting k = in Eq. ([1]). Since ((Tj^)k is an odd func- 
tion of ky as shown in Eq. the pseudospin in Fig. [1] 
at {kx, ky) and that at (k^, —ky) point to different orien- 
tations with respect to (<Ty). In contrast, the pseudospin 
at (kx,ky) and that at {—kx,ky) point toward the same 
orientation. Thus, the pseudospin component perpen- 
dicular to the zigzag edge flips, while the pseudospin is 
invariant for the armchair edge. 

We have seen for the zigzag edge that the reflection 
is intravalley scattering and that the pseudospin com- 



ponent perpendicular to the edge flips. For the armchair 
edge, the reflection is intervalley scattering and the pseu- 
dospin is invariant. More details concerning the scatter- 
ing, for example, the relative phase between the incident 
and reflected waves and the edge states are difficult to 
obtain within the above argument. In subsequent sec- 
tions we will explore an effective Hamiltonian to obtain 
the standing wave and the edge states. 



III. DEFORMATION-INDUCED GAUGE FIELD 

The fact that the pseudospin flips at the zigzag edge 
leads us to consider a gauge field for the edge that couples 
with the pseudospin in a manner similar to that an elec- 
tromagnetic gauge field couples with the real spin. Here, 
we show the formulation, in which the effect of the edge is 
included into the Hamiltonian as a deformation-induced 
gauge fieldi = 

To begin with, we consider a change of the nearest- 
neighbor hopping integral from the average value, —70, 
as — 7o+<57o,a(i"), where a {— 1, 2, 3) denotes the direction 
of a bond parallel to Ra in the inset of Fig. [T] The 
deviation 570, a (r) represents a lattice deformation in a 
graphene sheet. The low energy effective-mass equation 
for deformed graphene is written as 



H{v) 



*K(r) 
'J'K'(r) 
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/^*K(r) 
^*K'(r) 



(4) 



where ^'K(r) and vE'K'(r) are two-component wavefunc- 
tions that represent the electrons near the K and K' 
points, respectively. The Hamiltonian for deformed 
graphene is written as^i 



H{y) = i;f 



<T.(p + Aq(r)) 



<l>%r)ax 
<T'.(p-Aq(r))/ ' 



(5) 



where p = — i?iV is the momentum operator, a = 
{axjOy), and cr' = {—ax,Oy). A lattice deformation 
<^7o,a(r) enters the Hamiltonian through the deformation- 
induced gauge field A^(r) = {A%{y), A^(r)), where A'i(r) 
is expressed by a linear combination of (57o,a(r) as^i"— 



«F^2(r) = '57o,i(r) - i ((S7o,2(r) + 57o,3(r)) 



vyAI{v) 



V3 



(6) 



('57o,2(r) - (57o,3(r)) ■ 



The A'^(r) field causes intravalley scattering, while the 
perturbation that is relevant to intervalley scattering is 
given by a linear combination of A'^{r) and A'^{r) as2i 



<^'i(r) = (A2(r) + a^(r))e-2*'=^". 



(7) 



In Fig. [2ja), we consider cutting the C-C bonds lo- 
cated on the a;-axis at y = in order to introduce the 
zigzag edge in a flat graphene sheet. After cutting the 
bonds, the graphene sheet splits into two semi-inflnite 
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parts: y > and y < 0. The cutting is represented as 
i57o,i(i")ly=o = 7o, '57o,2(r) = and S'^oA^) = 0. From 
Eq. dH), the corresponding deformation- induced gauge 
field is then written as A'^{r) — {A'^{y),0), where A'^{y) 
is not vanishing only for the C-C bonds located on the 
X-axis at y = as A§(j/)|j,=o — (70 /^f)- Since ^§(y) is 
defined for the C-C bond, A'^{y) is meaningful when it 
is integrated from —^g to ^g, where is of the same or- 
der as the C-C bond length and will be taken to be zero 
at the end of calculation in the continuum limit. Note 
also that the vector direction of A'^(r) is perpendicular to 
that of the bond with a modified hopping integral. Since 
the zigzag edge is not the source of intervalley scattering, 
intervalley scattering can be ignored. Hereafter, we con- 
sider the electrons near the K point for the zigzag edge. 
Moreover, separation of variables can be employed, due 
to translational symmetry along the a;- axis. As a result, 
^K(r) and Px in Eq. ([5]) can be replaced with e*'^^^4'K(2/) 
and Px- The energy eigenequation can then be simplified 
as HYi{y)'^K{y) = E'^Kiy), where the Hamiltonian is 



HK{y) = VF [(Jx{px + A'^iy)) + CTyPy] ■ 



(8) 



This Hamiltonian is solved in Sec. IIVI The cutting 
which produces the Klein edges^li^^ is represented as 

<57o,i('') = 0' <^7o,2(r)|y=o = 70 and <57o,3(r)|y=o = 7o- 
From Eq. ([B]), the corresponding deformation-induced 
gauge field is then written as A'^(r) = (— Aj^(y), 0), where 
A'^{y) is the gauge field for the zigzag edge. Note that the 
direction of the A'i(r) field for the Klein edge is opposite 
that of the zigzag edge. 

The armchair edge can be introduced by cutting the 
bonds located on a; = ±6, as shown in Fig. [Ifb). By set- 
ting (57o4(r) = 0, (570,2 (r)U=-5 = 7o and 57o,3(i")U=5 = 
7o in Eq. ^ , the deformation-induced gauge field for the 
armchair edge is written as A'i(r) = (A§(a:), 0) with the 
limit of 5 — >■ 0. Due to translational symmetry along 
the 2/-axis, Py is replaced to Py in Eq. ([5|). Thus, the 
Hamiltonian is given by 



(Tx{Px + A'^ix)) + ayPy 

(jy^xYox 



(l>'^{x)cFx 

-CTxiPx - A'^ix)) + ayPy J ' 



with (^'^{x) = A%{x)e~'^'^^^^ . This Hamiltonian can be 
reduced further by means of a gauge symmetry in the 
following manner. Since A'^{x) does not depend on y, 
it can be represented in terms of a scalar function 1^9(2;), 
as A^{x) = dxfix). Using the gauge transformation: 



*K(a;) 



o{x) 



^Kix) and ^K'ix) 



A^{x) can be erased from the Hamiltonian for each val- 
ley. However, note that as a result of this gauge trans- 
formation, (/)'^{x) must be changed into e^*''''-^ To 
minimize notation, let us use to denote this gauge 

transformed field, so that ^^(x) = A'i{x)e^^^'^^'''>-''^''\ 
The Hamiltonian for the armchair edge is then written 



H{x) = Vf 



'JxPx+'JyPy 4>'^ix)ax 



<i}Hx)*ax 



(9) 



(a) zigzag edge 




(b) armchair edge 



X 








FIG. 2: (a) The bonds on the dotted line at 1/ = are cut 
to introduce the zigzag edge (Klein edge). The cutting is 
represented as a deformation- induced gauge field A''(r) = 
{A%(y),Q). (b) The deformation- induced gauge field for the 
armchair edge is given by A''(r) = {A'^(x),Q). 



This will be solved in Sec. [V] Note that by introducing 
Tq, (a — 1, 2, 3) matrices defined by 



/ 
/ 



T2 



-il 
il 



T-3 



/ 

-/ 



(10) 



the unperturbed Hamiltonian is represented in a compact 
fashion as -ffo(r) = vf{t3<7xPx + T^ayPy), where tq is a 
4x4 identity matrix. 

In Eqs. dH) and ([7]), we assume |(57o,a(r)| ^ 70, and 
ignore the higher order term of (57o,a(r)- As a result of 
this simplification, the relationship between A'^(r) and 
'^70, a (r) may deviate from Eq. ^ when |i57o,a(r)| ~ 7o- 
However, note that the direction and not the strength of 
the Ai(r) field can be determined by Eq. (jH]), even for the 
case where |(57o.a(r)| ~ 70. Consideration of this point is 
given in Appendix 1X1 
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IV. ZIGZAG EDGE 

The scattering problem for the zigzag edge is solved 
in this section. Standing wave solutions are constructed 
in Sec. IIV Al and the properties of the solutions are ex- 
amined in detail. Localized edge states are constructed 
in Sec. IIVBI The behavior of the standing wave in the 
presence of a mass term and an external magnetic field 
is examined in Sees. IIV CI and IIVDI respectively. The 
local density of states near the zigzag edge is calculated 
analytically in Sec. IIV El 



A. Standing Waves 

To begin with, solutions are constructed for the case of 
A^i^y) = in Eq. ©. Let $(?/) be the eigenstate of the 
unperturbed Hamiltonian H^{y) = vp (axPx + o'yPy) . 
H^{y) satisfies (TxH^{—y)ax = H^{y); therefore, a gen- 
eral solution may be constructed from the basis function 
$(2/) to satisfy the constraint equation. 



$(-2/) = e-'s^x^y), 



(11) 



where is a real number phase. The phase g can not be 
an arbitrary value. The successive operation of Eq. (|lip 
on $(—?/) gives $(—(—?/)) = e~^'3cr^$(?/), and hence 
g should be or TT. Note that a set of functions sat- 
isfying Eq. (fTTj) is useful for construction of solutions 
in the case A'^{y) ^ 0, because Hxiy) also satisfies 
<^xH},ii—y)crx = Hk(i/). This constraint comes from the 
inversion symmetry of the gauge field with respect to 

From Eq. ([TT|), we have $6(2/) = e*f$A(-2/). Thus, 
$(2/) can be rewritten as 



(12) 



By substituting Eq. ([12]) into H^{y)<^{y) = E^{y), we 
obtain simultaneous differential equations: 



E d 

— $s(y) = +Px^s{y) + h—^aiy), 

VF ay 

d 



lip 



(13) 



ay 



where ^s{y) and $a(y) are defined as 

$^(j/) EEE e~'*$A(2/) + e+'*$A(-2/), 
^aiy) = e^'i^Aiy) - e+'i^Ai~y)- 



(14) 



For the case g = 0, Eq. ([T4| implies that ^s{y) is an 
even function [$5(2/) = ^s{—y)], while $0(2/) is an odd 
function [$a(y) — —^ai—y)]- Thus, they can be param- 
eterized as follows: 



<^s{y) = Scos{kyy), 
<Pa{y) = Asm{kyy), 



(15) 



where the parameters S and A can be determined from 
Eq. (ini). By substituting Eq. ^ into Eq. ([13]), we 
obtain the secular equation 



= 0. 



(16) 



The solution of this secular equation satisfies 
E^^ihvFf{kl + kl), 



A 



-s. 



(17) 



Let ^(k) be the polar angle between vector k and the 
kx-SLxis. Then, kx = fccos0(k) and ky — A:sin0(k) 
where k = |k|, and the second equation of Eq. p7)) 
can be rewritten as A ~ 5tan [6'(k)/2] for the eigen- 
state with positive energy E = hvpk. Assuming that 
5 = cos [6'(k)/2], we have A = sin [6'(k)/2]. Substituting 
these into Eq. (IT5t gives 



^s{y) = cos (^^y^^ cos{kyy), 
^aiy) = sin sin{kyy). 



(18) 



Then, Eq. ([HI) is substituted into Eq. (HH) with g = to 
give 




Similarly, for the case where g = tt, we have 



(19) 



'sin (kyy - 
sin (kyy + 



(20) 



The energies of the eigenstates $"(2/) and $'^(2/) are 
equal, and therefore a general solution can be expressed 
as a superposition of the degenerate eigenstates, as 



$/(2/) = sin(/)$"(2;)+cos(/)<i>^(2/) 

'smikyy-e{k)/2 + f)^ 
^smikyy + e{k)/2 + f)^ 



(21) 



where / is a real number. The value of / is determined 
as follows. 

The Hamiltonian, H^iy) — H^{y) + v-paxA'^{y), is 
identical to the unperturbed Hamiltonian (y) for y to 
satisfy \y\ > ^g, so that $-^(2/) satisfies the eigenequation 
HK{y)^^{y) = E^f (y) for \y\ > ^g. We need to solve 
-ffK(2/)*K(2/) = E'i'Kiy) locally for \y\ < ^g. By parame- 
terizing the eigenstate of Hxiy) as 5'k(2/) — N{y)^^ {y), 
we obtain the constraint equation for N{y) and $-^(2/) as 

WyiPyNiy)] + axA2{y)N{y)} (y) = 0. (22) 
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To obtain Eq. ([22]) we must place ^'k(2/) = N{y)^f (y) 
to i?K(2/)*K(2/) = E^K{y), and use H°,{y)^f{y) = 
E^f{y). Here, we have assumed that the energy eigen- 
values of the standing wave ^K{y) and of are 
the same. This assumption is valid for the standing 
wave, because the energy eigenvalue is determined by 
the bulk Hamiltonian H^{y) and the energy does not 
change through elastic scattering. However, note that 
this assumption is not valid for the edge states, which 
HK{y)'^K{y) = E'^v.iy) must be solved directly (see 
Sec. lIVBl for more details). Now, Eq. ([^ is equivalent 
to the two successive equations: 



Al{y)N{y)-n'^^^ <i>^(y) = 0, 



(23) 



dy J 



The following two cases can be considered for this suc- 
cessive equation. One case is that the solution satisfies 

Al{y)N{y) + h^^^^Q, 

dy (l2/|<Cs)- (24) 

^ky) = o, 

The first (second) equation of Eq. ([^ ensures the second 
(first) equation of Eq. The other case is that the 

solution satisfies 



dy 



{\y\ < (25) 



AUy)Niy) 



The two conditions, Eqs. (|24p and ((25|) . correspond to the 
standing wave in the upper semi-infinite graphene plane 
for y > and that in the lower plane for ?; < in the 
limit of = 0, as shown in the following. 

For the case of Eq. (|24| , the first equation is integrated 
with respect to y, to obtain 



N{-Cg) = N{Qexp U Al{y)dy 



(26) 



Hence, when {l/U) /^^ A%y)dy » 0, N{(,g) is negli- 
gible compared with N{—^g), and therefore the stand- 
ing wave appears only for y < 0. In contrast, when 
{1/K) /^^ A'^{y)dy ^ 0, the standing wave appears only 
for y > 0. The other condition in Eq. ([M)) holds for 
the limit of ^ by setting / = -6'(k)/2 in Eq. ([2T|). 
because 



hm $r ''(''^/^(y) = 0. 



(27) 



This condition leads to ^k,b(0) — 0, which represents 
the boundary conditions for the zigzag and Klein edges 
shown in Fig. [^fa). Thus, Eq. (1^^ covers two situations, 
depending on the direction of the gauge field; A'^ (y) ^ 



or Al{y) < 0. That is, when A%{y) > 0, Eq. ([24]) cor- 
responds to the upper semi-infinite graphene plane with 
the zigzag edge, while when A'^{y) <C 0, Eq. (fM)) corre- 
sponds to the lower semi-infinite graphene plane with the 
Klein edge. Similarly, when A%{y) ^ 0, Eq. ([25]) corre- 
sponds to the lower semi-infinite graphene plane with the 
zigzag edge, while when A%{y) <C 0, Eq. (^5]) corresponds 
to the upper semi-infinite graphene plane with the Klein 
edge. 

From Eq. (pS)) . it follows that the gauge field for the 



edge should be large. 



{llh)fj, A%{y)dy 



> 1. In 



Ref. |29|, the following was obtained analytically 



1 



«9 



A2(y)dy = -ln(l-c) 



(28) 



where c is the parameter that specifies the deformation 
as (57o,i(i')|iy=o = c7o [see Fig. [2]. The right-hand side 
gives logarithmic singularities for c = 1 and c = —oo. 
The limit c — 1 corresponds to the zigzag edge, while the 
limit c — > — cxD represents the Klein edge. Note that when 
c — )• — oo, the electron is unable to have a finite amplitude 
on the A and B atoms located at y = 0, which effectively 
represents the Klein edge. Because of the singularity, 
N{ij) that satisfies Eq. ((26|) is similar to the step function; 
7V(y) = TV ^ for y < 0, and otherwise N{y) = 0. 

Now, by setting / = -6'(k)/2 in Eq. dUl), the standing 
wave in the conduction band is expressed as 



*K,k(r) = 



^N{y) 



/sin(fc,y-0(k))\ 
V sin(fcj,y) 



(29) 



where the plane wave parallel to the edge with the 
length Lj. is included. The standing wave in the va- 
lence band is obtained by using the particle-hole sym- 
metry of the Hamiltonian, <JzHK{y)'^z = ^H^iy), as 



n,k(y) 



^^*K.k(y): 



*K,k(r) 



^Niy) 



(sm{kyy-9{k))\ 
V -sin(fcj^y) J 



(30) 



Here, we consider the pseudospin of the standing wave. 
The pseudospin for an eigenstate 4'(y) is defined by the 
expected value of the Pauli matrices as (at) = J ai{y)dy 
{i = X, y, z), where cri(y) is a pseudospin density defined 
by cri(y) = '^Hy)'^i'^iy)- Note that the y-component 
of the pseudospin is proportional to the imaginary part 
of the Bloch function, such as (Jy{y) oc Im[\l/^^'B]. The 
Bloch function of the standing wave is real, so that the y- 
component of the pseudospin for the standing wave van- 
ishes, that is, ((Ty) = 0. Note also that (ay) = means 
that the current normal to the zigzag edge vanishes. 
It is interesting to note that (ctj,) — holds whenever 
^K,k,A(y) and 5'K,k,B(y) can be taken as real numbers. 
This indicates that the result (ay) = is not sensitive to 
the value of /, but depends only on the fact that does 
not have a relative phase between the two components. 
The condition of Eq. (I?71) means that the pseudospin den- 
sity is locally polarized into the positive z-axis near the 
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zigzag edge, that is, az{0) > and (JxiO) = cry{0) = 0. 
Actually, by substituting y ~ into Eq. (j^ . the stand- 
ing wave near the zigzag edge has amplitude only at 
A-atoms. This polarization of the pseudospin is con- 
sistent with the fact that the gauge field ^§(y) has a 
non- vanishing deformation-induced magnetic field, 



(31) 



at the zigzag edge. The presence of the B^{y) field at 
the zigzag edge causes local polarization of the stand- 
ing wave pseudospin near the zigzag edge, similar to the 
polarization of a real spin by a magnetic field. We will 
show in Sec. IIV El that this polarization of the pseudospin 
causes anomalous behavior to appear in the local density 
of states (LDOS) near the zigzag edge. 

A zigzag nanoribbon is given by introducing another 
zigzag edge at y ~ — L, in addition to the zigzag edge 
at y = 0. Suppose that the edge atoms at y = —L are 
B-atoms, which imposes the boundary condition on the 
wave function aty = —L as limj^=_L ^'k k Ai^) — ^- This 
leads to the constraint equation for (kx, ky), 



ky L 



-0(k) 



(32) 



where n is an integer. It is noted that this equation 
reproduces 



~kx tan(/ty^). 



(33) 



which was obtained by Brey and Fertig in Ref. [lo [the 
negative sign in front of kx is a matter of notation] . Note 
that n should be a nonzero integer, because the equation 
does not possess a solution when n — Q. For the case 
where the edge at y = — L is the Klein edge, the bound- 
ary condition on the wave function at y = —L becomes 
limy=_L 3(r) = 0. This leads to kyL — wk, where 
n is a positive integer. 



B. Edge States 

In addition to the standing wave derived in the previ- 
ous subsection, Hxiy) possesses localized edge states32. 
Here, we show how to construct the edge states. 

The following observation is useful in order to obtain 
the edge states. Instead of Eq. (|15l) . we assume 



$,(y) = 5cosh(2//0, 
$a(y) = Asinh(y/0. 



(34) 



By substituting Eq. into Eq. the secular equa- 
tion is obtained: 



s 



= 0. 



The solution of this secular equation satisfies 

= {hvpf (ki - r') , 



A 



-s. 



(35) 



(36) 



By introducing the </> variable, which satisfies 

= -kxtanhcf), (37) 
we have E'^ /{hvp)^ — k^J cosh^ 0. For the case 
E kx 

(38) 



TiVF cosh (j) ' 

we have A/ S = tanh(0/2). By inserting this into Eq. ((34l 
and setting S = cosh(0/2), we obtain 



^s{v) = cosh ( I j cosh J , 



$a(y)=sinh( -Jsinh|^- 



(39) 



(40) 



By substituting Eq. into Eq. dH]) with g = 0, we 

have 

\ cosh ( I — 
Similarly, for the case g — tt, we have 
/sinhff + 1 

-<^'H-(^.: 



(41) 



The energies of ^^{y) and ^^(y) are equal; therefore, the 
basis function may be chosen as 



(42) 



$+(y)^$"(y) + <i>-(y) = e+T 

<f-iy)^^"{y)~^^y)^e-^ (e+I/2) • 

The functions ^~^{y) and $~(y) are exponentially in- 
creasing and decreasing functions of y, respectively. 
Thus, neither ^~^{y) nor $~(y) is a normalized wave 
function all over the space, y € (—00,00). However, note 
that $^(y) and ^^(y) can be normalizable wave func- 
tions for y < and y > 0, respectively. We also note 
that the pseudospin of ^~^{y) is given by (az) = tanh0, 
while that of $~(y) is (cr^) = — tanhc/). 

From the above observation, we parameterized the lo- 
calized eigenstate as 

yi /p+g{y)\ 
^Kiy)^Ne-^(l_^^^,y (43) 

where A'^ is a normalization constant, and the modula- 
tion of the pseudospin is represented by a function g{y). 
Substituting Eq. (|^ into Hyi{y)^Yi{y) = E^y^iv) gives 
simultaneous differential equations for y(y). 



Px + Al{y) + h^ fM+5(y)') ^Z, 
dy \£, J VF 



,+2s(iy) 



(44) 
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By summing and subtracting both sides of Eq. (|44|) . the 
energy eigenequation can be rewritten as 



Px 



dy 



E 



cosh(25(y)), 



The sohition of the second equation is given by 



(45) 



9{y) 



[-ismh-i 


( hvF\ 








I hvp\ 


l^ + isinh-i 






kW) 



(y < 0), 
{y > 0). 



(46) 



The sign of g{y) changes across the zigzag edge, and this 
sign change indicates that the z-component of the pseu- 
dospin flips at the edge. The flip is induced by the gauge 
field A'^{y). To represent this, we integrate the first equa- 
tion of Eq. (pSj) from y — — to ^g, and acquire 



dgjy) 

dy 



dy 



Al{y)dy. 



(47) 



We have neglected other terms, because they are propor- 
tional to and become zero in the limit of — 0. By 
substituting Eq. (gH]) into Eq. (gT]), we find 



-sinh^if^ 
\^E 

lience, Eq. (gSj) becomes 



Al{y)dy. 



(48) 



9{y) = ' 



A''Jy)dy 



A'^{y)dy 



(y < 0), 
{y > 0). 



(49) 



Having described the wave function of the localized state, 
let us now calculate E and ^. To this end, we use the 
first equation of Eq. (|45|) for \y\ > and obtain 



E 

VF 



Px 



cosh 




(50) 



A''Ay)dy 



Moreover, using Eq. p8|) . we find 

r«9 




A^^{y)dy 



(51) 



In addition to this localized state, there is another local- 
ized state for the same with the same ^, but with the 
opposite sign of E. This results from the particle-hole 
symmetry of the Hamiltonian, and the wave function is 
given by a.^^K{y)- 



In the following, we will show that the solutions can 
reproduce all the properties of the edge states known in 
the tight-binding lattice (TB) model,— such as the 
asymmetric energy band structure with respect to the 
K (K') point, the fiat energy band, and the pseudospin 
structure. 

The asymmetric energy band structure with respect to 
the K (K') point originates from the normalization con- 
dition of the wave function, which requires that ^ should 
be positive. This requirement restricts the value of kx 
in Eq. (fSTI) . When A'^{y) is positive, Eq. (fSTI) indicates 
that the localized states appear only at fc^, < around 
the K point. This is the reason why the localized states 
appear in the energy spectrum only at one side around 
the K point. A similar argument can be used for the K' 
point, which concludes that the localized state appears 
at fca; > around the K' point. The Hamiltonian around 
the K' point is expressed by 



Hk'{v)=vpct' -{v^K^iv)). 



(52) 



Therefore, we obtain different signs in front of A'^ir) in 
H]^{y) and Eq. ([5^. which causes the negative sign in 
front of the right-hand side of Eq. ([5T|) to disappear for 
the K' point. Thus, when A'^{y) is negative (Klein edges), 
edge states appear on the opposite side; fc^; > around 
the K point and k^ <Q around the K' point. Calculations 
on the TB model with Klein edges also agree with the 
results obtained here. 

A singularity of the gauge field, |^§(y)| — > oo, 
is the origin of the flat energy dispersion and the 
pseudospin polarization of the edge states. When 
{l/K) /f^ Axiy)dy oo, E in Eq. (|50|) becomes zero. 
The zero energy eigenvalue between the K and K' points 
in the band structure corresponds to the flat energy band 
of the edge state^^i Moreover, from Eq. g{y) oo 
for y < and g(y) — >■ — oo for y > are obtained. 
In this case, the localized state is a pseudospin-up state 
\I'K(r) oc *(1,0) for y < and a pseudospin-down state 
^'K(r) oc *(0, 1) for y > 0. Hence, a singular gauge field 
at the zigzag edge causes polarization of the pseudospin 
of the localized states. Polarization of the pseudospin 
means that the wave function has amplitude only at the 
A (or B) atom, so that this result agrees with the re- 
sult from the TB model for the edge state. Comparing 
Eqs. jSI]) and ^ with Eqs. §7} and dSH]), the relation 
between the variable cj) and the field (y) is observed as 

cl>=il/n)j\ Al{y)dy. 

Here, we note that N in Eq. is not a function of y, 
but is a constant for y e (—00,00). Therefore, the edge 
states appear on both sides of the zigzag edge, y > and 
y < 0, while the standing waves appear on only one side 
of the edge under the hmit {l/K) jf^ A'^{y)dy ±00. 
With this hmit, the edge states can be confined to one 
side of the edge, because the energy of the locahzed 
state becomes E = 0, and therefore, the superposi- 
tion of an edge state, ^'k(2/), and its electron- hole pair 
state, (Jz'^Kiy), is a solution. It is easy to see that 
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^k(j/) + o'z^Kiy) has amplitude only for y < 0, while 
^k(2/) — f^z^Kiy) has amplitude only for y > 0. The 
wave function of the edge state for ?; < is then given by 



*K,fc.<o(r) = 



(53) 



where the normalization constant has been fixed, -y^ 2 1 /c^; | , 
by assuming that the system is a semi-infinite graphene 
plane. We note that the mass term, maz , is proportional 
to the particle-hole symmetry operator, az- Thus, the 
mass term automatically restricts the region where the 
edge states can appear {y > or y < 0), and this is 
shown in Appendix [Bl 

Finally, we consider the edge states in nanoribbons. 
Note first that the exact localization length for the case 
of a zigzag nanoribbon with width L satisfies 



,'L 

-krr tanh — 



o -k^L = (A;^^)atanh 



(54) 



which is obtained by analytical continuation ky — i/^ 
for Eq. (|33)) . Comparing this equation with Eq. ((5T|) 
shows that the large value of the gauge field in Eq. ([5T|) 
corresponds to the case oi L/^ ^ 1 in Eq. This 
is consistent with having solved the Hamiltonian locally 
near the edge, in which it was implicitly assumed that the 
condition L/^ ^ 1 is satisfied. Except the edge states 
whose localization length is in the order of L, Eqs. ([?T|) 
and ([M)) give almost identical values of ^ ~ which 
justifies the description using the gauge field. Note also 
that the condition L/^ ^ 1 also represents the condition 
kxL ^ —1, which is clear from the second equation in 
Eq. (|54l) . To solve the Hamiltonian for the edge states 
with ^ = 0{L), Eq. (gS]) must be solved globally, for 
example, on a circle, which is a challenging issue. 



Thus, by changing the basis state from j^) and j^) 

into |*K,k.A(B)) [= (l/y2)(|vE'^ k)± l*K,k))]. the Hamil- 
tonian is represented as 



( (*K,k,A|i?|?|*K,k,A) (*K,k,A|i/|?|* 
U^'K,k3l-ffKl*K,k,A) (*K,k,B|i/KN* 



K,k,B, 
'K,k,B' 



m hvpk 
hvpk —m 



Here, the angle (pk is defined as 
cos (pk 



m 



sm (pk 



fiv-pk 
Ek 



(57) 



(58) 



where Ek = \J vrP- -I- (hvph)-^. The normalized eigenvec- 
tors of the matrix in Eq. ([57)) are then 



(59) 



(^^^ <i^k I a,nd 
\sm ^ 



for the Ek and —Ek eigenvalues, respectively. 
^'K,k,A(2/) = V2N{y)sm{kyy - 9{k)) and 
^K.k.B(y) = V^N [y) sin{kyy) are obtained from 
Eqs. and ([50)1 : therefore, the standing wave near 
the zigzag edge is given by 



*K:k(2/) = ^N[y) 



's(%) sin(kyy-e{\<i))' 
sin sin {kyy) 

^ . sin ( ^ ) sin(fc,,j/ - 6'(k)) 
cos I ^ I sm [kyy) 



(60) 



where we have omitted to write the plane wave parallel 
to the edge. The factors cos (0^/2) and sin ((/'fe/2) appear 
in a manner similar to the eigenvalue problem of the spin 
magnetic moment in a magnetic field. 



Mass Term 



D. External Magnetic Field 



Let us reconsider the scattering problem for the case 
where the Hamiltonian includes a mass term. The total 
Hamiltonian is given by H^{y) = H}<i{y) + maz, where 
the mass, m, is a constant over the space y G (— oo, -l-oo). 
The solutions of H^{y) can be constructed from the so- 
lutions of HK{y) as follows. 

For HK{y), the standing wave solutions, Eqs. (1^ and 
((30| . satisfy 



^fK(2/)*K,k(y) = 7i^^Ffc*K,k(y), 



(55) 



For the mass term, because "^K.wiy) — '^z^'Kk(y)' 
obtain 



w^z*K,k(y) = "^*K,k(2/): 

maz'^Wiy) = m'^kkiy). 



(56) 



In this subsection, solutions are constructed for a mag- 
netic field applied perpendicular to the graphene plane^i 
A magnetic field B can be represented by the electromag- 
netic gauge field as A(j/) — (_By,0). This gauge field is 
included in the Hamiltonian by substituting the momen- 
tum operator p with p — eA. For the case ^^(y) = 0, 
the eigenequation becomes 



VF [(JxiPx - eBy) + ayPy] $(r) = E<^{r) 



(61) 



and the solutions are given by the Landau states, which 
are specified by an integer n and a center coordinate Y 
as 



$^^(r) = Cnvei^e 



^ l^sgn(n)y2Hi^|«|-l((^/-n/0^ (62) 

-H\n\iiy-Y)/l) J' 
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where CnY is a normalization constant, I = y^h/eB, 
and Hn{x) is a Hermite polynomial defined by Hn{x) = 
(-l)"e"= {d/dxYe-'-' [n > 0). The energy eigenvalue of 
$^y(r) is given by En = sgn(n)-\/2|n|?iwF/^- 

A method similar to that in Sec. IIV Al ls used to solve 
the scattering problem in the presence of a magnetic 
field. The energy eigenstate of Hk{v) is parameter- 
ized as ^'K(r) = N{y)^l;^{r). Substituting this into 
ifK(r)*K(r) = En^K{r), and using ffO(r)$LL(i.) ^ 
_E„$^y(r), we obtain the constraint equation for N{y) 
and $^L(i.), 

{<Jy[PyN{y)] + a^A2{y)N{y)} (j.) ^ q. (63) 

Two cases can be considered as a solution for this succes- 
sive equation [see Eqs. (fM)) and (f25|) ]. Here, we choose 
the case where 



A2{y)N{y) + h^ = 0, 



dy 



\y\ < 



(64) 



From Eq. (|62p . the second equation leads to 

C„ye~5(^)'iJ-|„|(_2y/0-0, (65) 

with the limit — >■ (y — >■ 0). The number of zeros 
of Hn is \n\, so that there are \n\ solutions of Eq. 
which are denoted as l^i (i = 0, • • • , |n|). The solutions 
can then be written as 



K,nYi 



,(r) = iV(y)$LL (r). 



(66) 



Note that 5'K,ny(i") with a large value of Y that satisfies 

Y ^ I, can be an approximate solution, due to the expo- 
nential factor in Eq. ((65)) . The solution with a large value 
of Y represents the wave function in the bulk, and is not 
sensitive to the details of the edge. The solutions given 
in Eq. ([66)) concern the Landau states near the zigzag 
edge, and these are examined in the following. 

For the case that n is an odd integer, y = satisfies 
Eq. ([S5)) . because 77„(0) = 0. The wave function with 

Y = decays according to exp{—y^/2P), and the ampli- 
tude has a maximum at the zigzag edge. Note that the lo- 
calization length is in the order oi I {I ~ 25[nm]/y^3[T]), 
which is larger than the localization length of the edge 
state ^ = — /CjT^ where takes a value of the same order 
as the lattice constant. 

The lowest Landau level {n — 0) can not satisfy the 
condition of Eq. ()65l) because Hq{—2Y/1) ~ 1 and the 
amplitudes of B-atoms do not vanish at the edge. Thus, 
the lowest Landau level is absent for the K point. On the 
other hand, the lowest Landau level appears for the K' 
point. The Hamiltonian for the K' point is given by 

i7K'(r) = «F [~cj.^{p^ - eA^{y) - Al{y)) -f ayPy] . (67) 

For the case of A'^{y) — 0, H]^{y) and Hy^i{y) are related 
as Hy^i{y) — ayHK{r)(Jy, and therefore the solutions for 



the K' point are given by (Ty<i>J^y, (r). The constraint 
equation for the K' point is 

{ay [pyN{y)] + cy.,Al {y)N{y)} ay<^'^, (r) = 0, (68) 

which reduces to the condition p^{y) = 0. The solu- 
tion is then given by 

^K'.nY^{v)^N{y)cyy-^^^,{v), (69) 
where Y,- denotes the solution of the constraint equation, 

C„ye-5(^)'y2Hff|n|-i i-y/l) - 0. (70) 

This condition is satisfied for n = 0, so that the lowest 
Landau level appears for the K' point. There is no con- 
straint for the value of Y . For the case of the first Landau 
levels (n = ±1), the Landau level for the K point appears, 
while that for the K' point disappears. Therefore, near 
the zigzag edge, the lowest and first Landau levels are 
not symmetric with respect to the K and K' points. 

E. Local Density of States 

Several groups have conducted scanning tunneling 
spectroscopy (STS) measurements to determine the 
LDOS near the step edge of graphite J^"— A peak struc- 
ture in the LDOS due to the edge states has been exten- 
sively discussed by many authors. Here, we calculate the 
LDOS near the zigzag edge. We show that some charac- 
teristic features that originate from the pseudospin po- 
larization, the edge states, and the mass appear in the 
LDOS. 

Let us first review the LDOS for graphene without an 
edge. Assuming that electrons are non-interacting, the 
bulk LDOS is given by 



1 



\E\ 



27r [fiVYY' 



(71) 



where p{E) is proportional to which results from 
the Dirac cone spectrum. Note that the actual LDOS is 
given by gsgvPiE), where gs = 2 [g^ = 2) accounts for the 
spin (valley) degrees of freedom. Next, the LDOS near 
the zigzag edge is calculated using the solutions given in 
Eq. dini)- The LDOS has the form. 



Ps{E,y) 



1 



\E\ 



R{E,y), 



(72) 



27r (hvp)^ 
where R{E, y) is defined as 

R{E, y)^- r d^?*™ (y)tvi;™ (y). (73) 

By performing the integral with respect to the angle 9 in 
Eq. (|73)) . we obtain an analytical result for R{E^ y) as 

' F{k\y\) + ^G{k\y\) {E > 0), 

R{E,y)^l (74) 
^ F{k\y\)~—G{k\y\) (S < 0), 
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where fc is a function of E according to fc = 
\/ E"^ — m? / (Tivp), and the functions F and G are defined 
as 



F{k\y\) = l- 



Jo{2k\y\) + J2{2k\y\) 



^^^1^1^ ^ M2k\y\) - J2{2k\y\) ^ 



Here, J,y{x) is a Bessel function of order v 

0.7 



(75) 



zigzag edge "'^^'^ 

y 




-0.6 -0.4 •. -0.2 • 

E (eV) 




energy E = — m should be polarized into B-atoms by the 
factors in Eq. ([5^ . and the amplitude of A-atoms are 
strongly suppressed by the mass term. Therefore, elec- 
trons with energy E = — m can not approach the zigzag 
edge, and therefore the LDOS disappears. Secondly, the 
LDOS peak of the edge states appears only at £' = 0.2 
eV. This is a straightforward consequence of the edge 
state amplitude appearing only for A-atoms. The ab- 
sence of the LDOS at i? = — m and the presence of the 
peak &i E = m due to the edge states occurs at different 
sides of the band edge. To plot the LDOS of the edge 
states in Fig. [31 we have used 



Pe{.E,y) 



1 



25 



1 



27r2 {E - my + 5^ 4y2 



(76) 



where (5 is a phenomenological parameter that represents 
the energy uncertainty of the edge states, for which we 
assume 5 = 10 meV. This result has been derived in 
Ref. Is^ for the case of m = 0. Note that pe{E,y) de- 
creases as ~ y~^, which is a slowly decreasing function 
compared with the exponential decay wave function of 
the edge state. 



FIG. 3: Positional dependence of the LDOS structure for the 
cases of m = 0.2[eV]. The number located on each solid line 
represents the distance (corresponding to \y\ in the inset) from 
the zigzag edge. The LDOS a,t E — —m vanishes near the 
zigzag edge, which is emphasized by the dashed circle. The 
dashed line denotes the LDOS in the bulk which is defined 
by the LDOS at \y\ — > oo. A peak structure due to the edge 
states is plotted for comparison. Note that there are several 
intrinsic perturbations^ that can change the position of the 
peak. 

Because the case of m = has been considered else- 
where;^ we consider the case to ^ here. Eq. ([7T|) holds 
for \E\ > \m\. The bulk LDOS vanishes for the case 
\E\ < |?7i|, as shown by the dashed line in Fig. [31 Note 
that the LDOS disappears suddenly a,t E = zLm, and the 
bulk LDOS has a step like structure at E = ±|to|. The 
bulk LDOS is symmetric with respect to E' = 0, even for 
the case m 7^ 0. However, note that the LDOS near the 
edge is not symmetric for the case of m ^ 0, which is 
clear from the different signs in front of the function G 
in Eq. In Fig. H the LDOS are plotted at \y\ = 0, 

1, 2, and 3 [nm] for the case of m = 0.2 eV. Note that for 
the case to = —0.2 eV, the corresponding LDOS curve is 
given by interchanging the conduction and valence bands 
in Fig. H 

The asymmetry in the LDOS near the edge appears 
at the following points. First, a step structure appears 
only at E = 0.2 eV. At E = -0.2 eV, the LDOS van- 
ishes, and the step structure is absent, as indicated by 
the dashed circle in Fig. [3l The absence of the LDOS at 
E = —0.2 eV can be explained by the zigzag edge consist- 
ing of A-atoms makes the standing wave polarized into 
A-atoms near the zigzag edge. However, eigenstates with 



V. ARMCHAIR EDGE 

In this section, the scattering problem for the armchair 
edge is solved using a method similar to that used in 
Sec. IIVI The standing wave shows that the pseudospin 
does not change its direction through the reflection at 
the armchair edge. 



A. Standing Waves 

Solutions for the case of (^'^{x) = in Eq. ([9]) are con- 
structed first, and then used as the basis functions to 
construct the standing wave near the armchair edge. Let 
$(a;) represent the solution of the unperturbed Hamil- 
tonian, Ho{x) = VFirsaxPx + '''Q'^yPy)- The perturbed 
Hamiltonian satisfies H{—x) — tiH{x)ti, so that the 
functions ^{x) that satisfy the constraint equation 

ri$(-a;) = e"'s$(x), (g = 0, tt) (77) 

are useful for construction of solutions in the case of 
(/)'^(x) 7^ 0. From Eq. (|77p . we may write 



<S>k{x) 



^e'9$K(-a;) 

By using Eq. (|75)) . the energy eigenequation becomes 



(78) 



f E . \ . , . . d 



(79) 
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where <I>s(x) and ^a{x) are defined as 



(80) 



For the case 3 = 0, $s(a;) is an even function, while ^a{x) 
is an odd function. For the case of 5 = tt, $s(x) is an 
odd function, while $a(a^) is an even function. 
For the case g = 0, we can set 



^six) = cos{kxx)(f>, 
^a{x) = ism{kxx)(l). 



(81) 



Substituting these into Eq. ([75)1 . we obtain the secular 
equation: 



E 



a- ■ k 



0. 



(82) 



The solutions of this secular equation satisfy E = ±hvpk, 
and the eigenfunction in the conduction band is given by 



''K.k' 



which is defined as 



V2 



9(k) 



1 



(83) 



By substituting Eq. (|8Tj) into Eq. (l80l) . we obtain 
$K(a;) = e**^-"^./)^ k/2- Using Eq. it can be seen 

that 



Similarly, for the case ^ = tt, we have 



(84) 



(85) 



New basis functions are defined using Eqs. (|84l) and (|85j), 

as 







$K(2;)^i(<i>"(x)+<I'-(2;))=0^,k i^' 



(86) 



The eigenstate $^(a;) represents a free propagating state 
with momentum k near the K point, while $^ (a;) repre- 
sents a state with momentum k' = {—kx,ky) near the K' 
point. It is clear that these are eigenstates in the absence 
of the edge. In the presence of the armchair edge, neither 
$^(a:) nor $^ {x) is an eigenstate, but a true eigenstate 
is the standing wave that is given by a superposition be- 
tween ^^{x) and (x) as 



*(x) ^ c'^{x)<^>'^{x) + c^' {x)<^^' (x) 



(87) 



To find c^'^ (x), it is useful to rewrite the total Hamil- 
tonian as 



H{x) = Ho{x) + VF<Tx [<P'^{x)ti + 0?(x)r2] , 



where (j)'^{x) is expressed in terms of real and imaginary 
parts, as (j)'^{x) = (j)'^{x) — i(f)1{x). In Sec. HID we have 
shown that 01 (x) = A;^(x)e2*[v'(^)-'='=-^l , where Al{x) = 
dx^p{x). A'^{x) is an even function with respect to x; 
therefore, (/'(x) can be taken as an odd function, so that 
the field satisfies (t>'^{—x) = (f>'^{x)*. From this condition, 
it follows that (j)'^{x) is an even function, while (j)l{x) is 
an odd function. 

Next, we construct solutions for the case (t)'^{x) — 0. 
Let us define (pi{x) and ^P2{x) using a real function f{x) 
as 



ipi{x) 

IP2{X) 



/cosh/(x) sinh/(x) 
lsinh/(a;) cosh/(x) 



(89) 



Since $^(x) and $^ (x) are the solutions of Ho{x), we 
obtain the following equations for f(x) from H{x)ipi{x) = 
E(pt{x), 



Px -i(l)i{x)\ f cosh f{x) sinh/(a:) 



i4>t{x) -px J \sinh/(x) cosh/(x) 



0. (90) 



All four components of this matrix are reduced into the 
same differential equation: dxf{x) = —(j)i{x)/h. 4ii{x) is 
an odd function, so that we have /(^g) = f{~S,g) by using 
/^^ 4>i(x)dx = 0. Because f{x) = when 4>i(x) = 0, 
the constant of integration can be taken as zero. As a 
result, we have /(^g) = f{^£.g) ~ 0. Therefore, f{x) can 
take only a non zero value for \x\ < ^g, and the mixing 
between $^(a;) and $^ (x) is negligible in the bulk. 

Finally, we assume that the solution of the total Hamil- 
tonian of Eq. ([88l) has the form of 



*±(x) =iV(x)[(^i(x)T»<P2(x)]. (91) 
The constraint equation for N{x) is then given by 

. dN{x) 



4>r{x)N{x) ± n- 



dx 



: 0. 



(92) 



This constraint equation has the same form as Eqs. ([M)) 
and (j25p . Performing the integral for x from — to ^g 
in Eq. dH!]) gives 



exp 




(j>r{x)dx 



(93) 



As we have shown in Fig. [5J (j)r{x) is a negative large 
quantity. Thus, ^-|-(x) has an amplitude only for x > ^g, 
while ^'-(x) has an amplitude only for x < — Cg- ^1(2^) = 
$^(x) and ip2ix) = (x) for |x| > ^g; therefore, the 
standing wave near the armchair edge is written as 



*k(r) 



gikyy 



N{xWk.^ 



— ikx 



(94) 



Note that the Bloch functions for the K and K' points 
are the same, which indicates that the pseudospins of 
the incident and refiected waves are equal, as shown in 
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Fig.m Thus, the Berry's phase of the standing wave near 
the armchair edge is given by — tt, which is in contrast 
to the case of the zigzag edge.— The boundary condition 
for the armchair edge does not forbid an electronic state 
to cross the Dirac singularity point, and therefore the 
electron can pick up a nontrivial Berry's phase. 



reflected wave 




K 

incident wave 



FIG. 4: The armchair edge reflects the wave vector k = 
{kx,ky) of one valley into k' — {—kx,ky) of another valley, 
and the two wave functions of the different valleys form a 
standing wave. The pseudospin is unchanged by the armchair 
edge. Note that the pseudospin for states near the K' point is 
not parallel to the vector k', while the pseudospin for states 
near the K point is parallel to the vector k. 

To understand the behavior of the standing wave in 
more detail, the density of Tq, was examined. The den- 
sity for an eigenstate ^'(j/) is defined by the expected 
value of ) = '^^{x)Ta'^{x). It is then straight- 

forward to check from Eq. ([M)) that ti{x) (x ± sm{2kxx), 
T2{x) oc zLcos{2kxx), and T3{x) = 0. ti{x) vanishes near 
the armchair edge (at a; = 0), and T2(x) takes a maximum 
value at the edge. This behavior can be understood from 
Eq. (I55|) . in which n couples with 4>r{x). Since (l)r{x) is 
singular at a; = 0, ti{x) can not have a non-zero value at 
X = 0. The result T3{x) = indicates that time-reversal 
symmetry is preserved. 



B. External Magnetic Field 

Let us examine the Landau states near the armchair 
edge. The electromagnetic gauge field A = (0, —Bx) for 
an external magnetic field B is included in the Hamilto- 
nian H{r) of Eq. ([5]) by the substitution p — >■ p — eA. 
The Hamiltonian satisfies H{—r) = TiaxH{r)axTi, and 
therefore the solution can be written as 



*K(r) 
?a,*K(-r) 



(95) 



Let "I>K(r) be the solution for the case (j)'^{x) = 0. Then 
<f'K(r) satisfies the following energy eigenequation: 

VF [<JxPx + cTyiPy + eBx)] $K(r) = E<^K{r). (96) 



The solutions are the Landau states specified by integer 
n and a center coordinate X as [see Eq. ([S^ ] 



sgn(7i) 




H_i((x-X)/0 
i{x-X)/l) 



(97) 



Applying the parity transformation r — — r to ^l^nxi^)^ 
we obtain 



LL 

K.nX 



i-r) = (-1)"+V,<i> 



K,n-X 



(r). 



(98) 



The matrix cr^ on the right-hand side can be understood 
by applying the parity transformation r — > — r to this 
energy eigenequation: 

VF [cTxPx + CTyiPy + eBx)] $K(-r) = -£:<I>K(-r). (99) 

The negative sign in front of the right-hand side shows 
that the energy eigenvalue of $k(— r) is opposite to that 
of $K(r)- By substituting Eq. ([Ml) into Eq. (US]), we 
obtain 



^tnxir) 



(100) 



By repeating the same argument given in the previous 
subsection, the following standing wave solutions are ob- 
tained: 



LL 
nX± 



(x) = N{x) 



(101) 



There are no constraints for the value oi X. It is then a 
straightforward calculation to check that ti{x) vanishes 
at the armchair edge. 



VI. DISCUSSION AND SUMMARY 

A realistic graphene edge may be a mixture of zigzag 
and armchair edges.J^^--^ The construction of the stand- 
ing wave near the general edge is one of the interest- 
ing applications for our framework. We believe that the 
Hamiltonian in Eq. ([5]) can describe the low-energy elec- 
trons in a graphene plane with a general edge. How- 
ever, note that this issue is related to the coherence 
length of the standing wave. In the present paper, we 
have not considered perturbations that break coherence, 
such as electron-phonon interaction. Interestingly, the 
electron-phonon interaction can also be represented as 
a deformation-induced gauge fieldi^i^ Thus, the gauge 
field description for the graphene edge may be useful 
when we consider such issues. 

The effective-mass model of Eq. ([5]) is equivalent to a 
chiral gauge theory for graphene that has been proposed 
by Jackiw and Pi.-^ Indeed, by applying ax to \1/K'(r) in 
Eq.(|4]), the Hamiltonian in Eq. ([5]) may be rewritten as 



H' = VF 



a- ■ {p + A'i{r)) 
0q(r)* 



0<i(r) 
-*T.(p-A<i(r))^ 



(102) 
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which is the electronic Hamiltonian of the chiral gauge 
theory. They have investigated zero-mode solutions of 
the Hamiltonian with a topological vortex for A'i(r) on 
the background of Kekule distortion for (p"^ (r) , in the con- 
text of fractionalization of quantum numberi^""— Our 
trial is then to study the graphene edge as a chiral gauge 
theory, although our results in this paper do not clar- 
ify fully the topological features of the graphene edge. 
It is interesting to note that one may find an advan- 
tage of a chiral gauge theory when we consider the real 
spins of the electrons. For example, the magnetism of the 
edge states may be understood as a parity anomaly phe- 
nomenoni^i^ The various field-theoretical techniques 
may be utilized to explore the electronic properties near 
the edge. Note also that the perturbation which mixes 
the electrons in the two valleys has been examined in the 
studies on the topological defect in graphene . ^^'^^ 

We have taken into account the edge as a part of the 
Hamiltonian. This strategy stems from the tight-binding 
lattice model, in which the edge is automatically included 
as a part of the Hamiltonian. A similar concept is found 
in the article by Berry et al}^ , in which the authors mod- 
eled the edge using a mass term, V{r)az- They consid- 
ered that a singularity of the mass V{r) — )• oo outside 
of the edge is necessary, in order to uniquely specify the 
pseudospin. We have observed a similar situation for the 
deformation-induced gauge field for the edge, that is, the 
field is singular at the edge. It is also interesting to note 
that (j)'^{r) in Eq. (|102l) corresponds to the mass of a Dirac 
fermion, and that the armchair edge is a singular point 
as for the mass. 

In summary, we have proposed a framework in which 
the edge is represented as the deformation-induced gauge 
field. We have used the framework to investigate the 
standing waves and edge states in the presence of a mass 
term and a magnetic field. The description of the edge 
using the deformation-induced gauge field is one attempt 
to better understand the edge. If we can describe the 
variety of edge structures as different configurations of a 
single gauge field, it provides a basis to further explore 
the properties near the edge. 
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Appendix A: Rotation of Pseudospin 

The configurations of the pseudospin field for three 
equivalent corners of the graphene BZ are not the same, 
as shown in Fig. [T] Consideration of this pseudospin be- 
havior is given in this Appendix. 



The tight-binding Hamiltonian can be written as2i 

where /a(k) ee e*'^" (a = 1,2,3). Note that /a(k) sat- 
isfies 

/a(k + nbi +mb2) = /,(k)e-^^("+'"\ (A2) 

where n and m are integers. Hence, the representations 
of H{k), H{k+hi), and i7(k-|-bi -|-b2) are different from 
each other, and are related via ff(k-l-bi) = MH{'k)AI~^ 
and H{k + hi+ bz) = M-^H{k)M, where 

/ +i2^/3 \ / 27r \ 

M=[\ gJvaj =oxp[^-a.y (A3) 

For the solution of H(k), we have the correspond- 
ing solution of the effective Hamiltonian at k -I- bi as 
^'bi = M"^. M is a rotational matrix for the pseudospin 
around the z-axis, so that the pseudospin of 4* and that 
of M'^ are related by rotation around the z-axis by an 
angle of 27r/3. This explains why the configurations of 
the pseudospin field around the three equivalent K (K') 
points are different from each other, as shown in Fig. [TJ 
Next, we consider the effective Hamiltonians for three 
equivalent K (K') points. By expanding /a(k) around 
the wave vector of the K point, kp — (47r/3a,0), we ob- 
tain faikp +k) ^fa{kF)+ifa{kF)k-Ra-\ . Usiug 

/i(kF) = 1, /2(kF) = e~'^, and /3(kF) = e+'^, we 
have H{kF -f k) = vpa • p + • • • , where p = ?ik and vp — 

37oacc/2?i. Then H{kp +hi + k) = Mvptr ■ pM'^ -) , 

and -ff(kF + bi + ba + k) = M'^vpcr ■ pM + ■■ ■ are 
obtained. The same argument can be applied to the 
K' points. For the K' point at — kF, we obtain the ef- 
fective Hamiltonian iJ(— kF + k) = vpcr' • p + • • • . It 
is then straightforward to obtain _ff (— kF + bi + k) = 
Mvpcr' ■ pM~i + • • • , and H{~kp -f- bi + b2 k) = 
M^^vpcr' ■ pM -I- • • • . This difference in the representa- 
tions of the effective Hamiltonians does not cause a prob- 
lem, because a coordinate transformation can be used to 
eliminate the M matrix from one effective Hamiltonian 
(see also Appendices in Ref. 111)4^ Here, we imply the 
coordinate transformation as the rotation of the x and y- 
axes by ±27r/3. A coordinate transformation cannot al- 
ter the physics, and therefore the physical result derived 
from the effective Hamiltonians are the same. Rather, 
by using the change of the effective Hamiltonians under 
a translation given by the reciprocal lattice vectors, a 
constraint for the form of the effective Hamiltonians can 
be obtained. For example, the deformation Hamiltonian, 
cr • A'^(r), should transform in the same way as cr • p. 
Therefore, we must have H{kp+k) = vpcr-{p+A'^) + - ■ ■ , 
H{kp +hi +k) = MvpCT ■ (p -I- A^)M-i + • • • , and 
ff(kF-f bi-hb2-Fk) ^ M-^vpcr ■ {p + A'^)M + ■ ■ ■ . Oth- 
erwise, there would be three physically distinct effective 
Hamiltonians for the same K point. The deformation- 
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induced gauge field Ai(r) satisfies this constraint, be- 
cause 



t,F(AS(r) -zA^(r)) = ^<57o,a(r)/.(kF) 



(A4) 



Note that this equation is equivalent to Eq. ©. The 
phase factor of e^*^'^/'^ appears when we change kp to 
kp + bi and to kp + bi + b2, due to the factor of /a(kp) 
on the right-hand side. A notable feature is that the con- 
straint must be satisfied for a strong lattice deformation 
that corresponds to a large value of A'i(r). Therefore, the 
direction of the gauge field does not change, although the 
values of 70, a (r) are renormalized for a strong deforma- 
tion. 



Appendix B: Edge states and Mass 

The edge states in the presence of a mass term is of 
interesting, because the magnetism of the edge states is 
related to the generation of a local spin-dependent mass 
term due to the coulombic interaction;^ Here, we show 
how to obtain the edge states in the presence of a uniform 
mass term. 

By substituting Eq. gS^i into H^{y)-^K{y) = £'*k(?/), 
we obtain instead of Eq. (j45|) 

. dg{y) 



Px 



dy 



Dcoshi2giy) + f), 



(BI) 



where the variables D and / are respectively defined as 
D EE ± — y^E^ - m2 and tanh(/) = (B2) 

Dp E 

The solution of the second equation in (IBip is 

n 



My) + f 



sinh 



sinh 



h 



(y<o), 
{y > 0). 



(B3) 



The first equation in Eq. (|B1|) is integrated with respect 
to y from —^g to ^g. Considering the limit ^g 0, only 
singular functions of ^^(y) and g{y) at y — can survive 
after the integration, so that we obtain 

r«9 



sinh 



h 
Id 



AUy)dy. 



(B4) 



Using Eqs. (jB3|) and (|B4p . we see from the first equation 
in (Ell that 



D = 



Px 



cosh 



(B5) 



holds except very close to the edge. From Eqs. (jB4p and 
(|B5p . we see that ^ is given by 



— — —kr tanh 




Al{y)dy 



(B6) 



Note that ^ in the presence of a mass term is identical 
to f in the absence of the mass given in Eq. ([5T|) . The 
mass term would affect ^, but this is not the case. From 
Eqs. (|B2[) and (jBSP we obtain the energy eigenvalue 



E = ±^Jm^ + {vpOy. 



(B7) 



When /^^ A°^{y)dy -> 00, we obtain D = and E = 
±|to|. 

According to the definition of / in Eq. (IB2p . the sign 
of / depends on the signs of both m and E. Let us 
first consider the case of i? < 0, by which we have / = 
sign(m)|/|. Using this expression for / in Eq. (|B3I) . we 
obtain 



5(2/) 



^ A?(y)d2/-isign(m)|/| (y < 0), 

1 1 

AUy)dy - i:sign{m)\f\ (y > 0). 



(B8) 



To determine |/|, we substitute Eq. (IB7| into Eq. (|B2l) . 
and considering that tanh(|/|) can be approximated as 

l-2e~2|/l for I/I > 1, we then have |/| « f^'^^ A'^{y)dy 

for A'i(y)dy > 0. Since /f^^ A'^{y)dy > for the 
B8I) becomes 



/-^ A%{y)dy 
zigzag edge, Eq. ( 




Al{y)dy (y<0), 
(y > 0), 



(B9) 



when TO < 0. The wave function of this eigenstate for 
?/ > 0, which has unpolarized pseudospin, is negligible 
due to the normalization. Thus, the localized state with 
energy _B = — |to| in the valence energy band can appear 
near the edge only for j/ < 0, and the wave function is 
given by ^K{y < 0) oc exp(— |y|/^)*(l, 0). Similarly, for 
the case of to > 0, we have 



9{y) 



iy < 0), 
A2iy')dy' (y>0). 



(BIO) 



The corresponding wavefunction is has the pseudospin 
down state, which appears only for y > near the edge. 
It is noted that the mass term automatically selects the 
region where the edge state can appear, y > or y < 0. 
This is reasonable, because we have used the particle-hole 
symmetry operator CTz to restrict the edge state only for 
?/ > or y < in Sec. IIVBI The particle-hole symmetry 
operator is nothing but the mass term. 
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